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Magnetic resonance imaging method 

[0001] The present invention relates to a magnetic resonance imaging method 
to produce successive magnetic resonance images. 

[0002] MR imaging has been established as a viable and unique imaging 
technique for cardiological examinations. Its utility is related to the ability to provide 
both precise morphology and assessment of global and regional ventricular function. 
For functional assessment of the heart, steady-state free*precession (SSFP) has 
become an indispensable sequence due to its superb blood-to-myocardium contrast 
(see e.g. Oppelt A. et al. FISP: a new fast MRI sequence. Electromedica 1986; 54: 
15-18.). Compared to spoiled steady-state sequences, SSFP also offers improved 
signal-to-noise ratio (SNR), and it allows for a shorter minimum repetition time (TR). 
However, the advantages of SSFP come at the expense of increased susceptibility 
to experimental non-idealities. As a result, it is only with recent advances in scanner 
hardware and sequence design that SSFP has become widely available. One 
prominent artifact in SSFP imaging is the so-called banding artifact. This results from 
unaccounted dephasing of the magnetization, such as from inhomogeneity of the 
main magnetic field. To reduce susceptibility to banding artifacts, it is critical to keep 
TR as short as possible, so that there is less time for unaccounted dephasing to 
occur. This, in turn, limits the use of echo-planar readout and, accordingly, 
acquisition efficiency of SSFP imaging is lower than in multiplanar image formation 
(see Mansfield P. Multiplanar image formation using NMR spin echoes. J Phys C 
1977; 10: L55-58.). 

[0003] SSFP can be applied in a real-time mode or in a triggered mode. In real- 
time mode, data are acquired in a continuous fashion. The lower acquisition 
efficiency of SSFP limits the achievable spatial and temporal resolutions. In the 
triggered mode, the idea is to "freeze" the effects of anatomical motion, which would 
reduce image quality. This is achieved using cardiac triggering - the coupling of 
image acquisition to a trigger signal - for cardiac studies, generally to an ECG signal. 
Using the assumption that most biological motion is rhythmic, triggering permits data 
acquired at different cardiac cycles but at the same point within the cycle to be 
combined to improve signal-to-noise ratio in the final image or to aggregate sufficient 
data to reconstruct a fully-sampled image. Since the images from multiple cardiac 
cycles are acquired at the same point within the cycle, the effects of motion are 
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minimized. However, residual motion does exist, due to cycle-to-cycle variation, and 
it is manifested as image artifacts. Therefore, it is desirable to reduce the overall 
scan duration in order to minimize such variability. A series of Images acquired at 
different points of the cardiac cycle can also be played back in cinematic loop 
fashion, called cine imaging, to provide physicians a way of viewing dynamic 
activities. Triggered acquisitions are generally performed In conjunction with a fast 
Imaging pulse sequence for Improved temporal resolution. The lower acquisition 
efficiency of SSFP mentioned previously leads to a longer overall scan duration, 
compared to echo-planar imaging sequences, for example. 

[0004] An Important application of SSFP Imaging is In obtaining volumetric 
time-resolved images of the heart. For example, multi-slice, multi-breathhold imaging 
based on SSFP is currently the method of choice for obtaining cine views of the 
heart. From these images, functional parameters can be determined with high 
accuracy and reproducibility (see e.g. W. Li et.al. MR assessment of left ventricular 
function: quantitative comparison of fast imaging employing steady-state acquisition 
(FIESTA) with fast gradient echo cine technique. J. Magn. Reson. Imaging, 2002; 16, 
p. 559-564). There are, however, limitations associated with multi-slice, multi- 
breathhold acquisitions. Inconsistencies in the multiple breath-holds may lead to 
misregistration of slices. Furthermore, the low bandwidth-time product of the 
excitation pulses used for short TR SSFP imaging results In large variations of 
excitation angles within the slice. This not only affects the signal-to-nolse and 
contrast-to-noise ratios, it also hampers the reconstruction of multi-planar reformats 
through the heart due to the slab boundary artifact. 

[0005] To overcome the aforementioned shortcomings, SSFP with true 3D 
imaging was previously developed (see B.A. Jung et. al. Single-breathhold 3D-true 
FISP cine cardiac imaging. Magn. Reson. Med. 2002; 48. p. 921-925). Cine 3D 
SSFP of the heart was shown to offer superior signal-to-noise and contrast-to-noise 
ratios compared to multislice sequences. However, in that work, compromises 
between Image quality, spatial and temporal resolutions had to be made to 
accommodate sufficient volume coverage Into a single breath-hold. Using an 
elliptical shutter in ky-kz space, cine volumetric data sets may be obtained at 
relatively high spatial resolution in a single breathhold. Still, temporal sampling does 
not satisfy the requirements for accurate left ventricular volume calculations. 
Generally, at least 10 temporal frequency harmonics are considered sufficient for 
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accurate volume curve calculation, requiring a sampling rate of at least 20xHR/60 
Hz, with HR being the actual heart rate. 

[0006] It is an object of the present invention to speed up SSFP dynamic 
imaging, particularly by a high acceleration factor. The acceleration factor is defined 
as the ratio between the acquisition time of a conventional scan with full sampling of 
the data to the acquisition time of an accelerated scan with reduced data sampling. 
[0007] This object is achieved by an MR imaging method with the features of 
claim 1. A practical use of this invention, as described in detail below, is in 
accelerating cardiac cine 3D imaging in order to obtain a sequence of cine images in 
a much shorter time or with a higher temporal resolution for a fixed acquisition 
duration. Another promising use of this invention is in accelerating free-breathing, 
untriggered (referred to as real-time hereafter) imaging to achieve higher spatial 
and/or temporal resolutions. 

[0008] The present invention is based on the knowledge that the high intrinsic 
signal-to-noise and contrast-to-noise ratios of SSFP imaging make it a good 
candidate for reduced data acquisition schemes. Moreover, sparse sampling of /c- 
space may allow shortening of scan time. At the same time, shorter scan time have 
the added benefit of reducing the overall power deposited in the subject, which is a 
concern for cardiac SSFP imaging especially at higher field strengths. 
[0009] However, accelerated SSFP is prone to image artifacts resulting firom 
instabilities of the gradient hardware system, such as eddy currents. The gist of this 
invention is to modify the profile order (also known as the phase-encode order) to 
reduce the said artifacts. These artifacts are often present in SSFP, but they 
generally remain at low level in practice. Therefore, linear profile order sampling in k- 
space prevails, particularly in 2D imaging. In the case of high acceleration factors, 
which are desirable in terms of reduced scan time or increased spatiotemporal 
resolutions, these artifacts are amplified drastically. In those cases, the resulting 
image artifacts overwhelm the image contents and must be suppressed. 
[0010] Further advantages of the invention are disclosed in the dependent 
claims and in the following description in which an exemplified embodiment of the 
invention is described with respect to the accompanying drawings. It shows 
[0011] Fig. 1 regular sampling pattern of magnetic resonance signals. 
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[0012] Fig. 2 time series of fully sampled images 1 and corresponding signal 
distribution of object intensities 3 along the dashed line 2 after inverse Fourier 
transform over time. 

[0013] Fig. 3 schematics of the k-t BLAST method, 

[0014] Fig. 4 Cine 3D SSFP: influence of sampling order on image artifacts, 
[0015] Fig. 5 Real-time 2D SSFP: influence of sampling order on image 
artifacts 

[0016] Fig. 6 the sampling pattern for cine 3D k-t BLAST with integrated 
training data, 

[0017] Fig. 7 a cine 3D single breath-hold data set obtained with 5x k-t 
BLAST in a volunteer, 

[0018] Fig. 8 a cine 3D single breath-hold data set obtained with 5x k-t 
BLAST in a cardiac patient, and 

[0019] Fig. 9 diagrammatically a magnetic resonance imaging system. 
[0020] The MR method according to the present invention is based on an 
approach for accelerated data acquisition in dynamic imaging, called k-t BLAST 
(Broad-use Linear Acquisition Speed-up Technique) and k-t SENSE (SENSitivity 
Encoding). The k-t BLAST and /c-f SENSE methods are described in an abstract byJ. 
Tsao et. al. (Prior-information-enhanced Dynamic Imaging using Single or Multiple 
Coils with k't BLAST and k-t SENSE. Proc Int Soc Magn Reson Med, Abstract 2369. 
2002). 

[0021] In the following, the k-t BLAST method is briefly reviewed, k-t SENSE 
is similar, but it additionally takes the complex sensitivity information of the receiver 
coils into account, and it is described in detail in the aforementioned document. 
[0022] Image content in a cine series of the heart exhibits a high degree of 

correlation in both space and time. This is demonstrated in Figures 1 and 2. Figure 1 
shows the sampling pattern in k-t space {k = wavenumber, t = time), indicating that k- 
space is fully sampled at every time point t The density and extent of sampling in /c- 
space are chosen according to the desired field of view and resolution, respectively. 
Figure 2 shows the corresponding fully-sampled, cine short-axis acquisition of the 
heart 1. Image information along a vertical profile 2 through the heart (Figure 2) is 
confined to compact regions in x-f space 3 (x = spatial position, f = temporal 
frequency), according to the dynamics of the corresponding image pixels. In the k-t 
BLAST method, acquisition efficiency is increased by applying sparse sampling to k-t 
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space according to a sheared grid pattern 5 (Figure 3c). This results in a tighter 
packing of signal replica in the reciprocal x-f space 6. Using low-resolution image 
information sampled for all cardiac phases 7 (Figure 3a), an estimate 8 of the 
expected signal distribution in x-f space is obtained, which is used to determine a 
reconstruction filter 9 (Figure 3b). The sparse sampling in the actual, high-resolution 
acquisition results in multiple replica 6 in x-f space, which in turn manifest as multiple 
fold-over in the actual image series 10 (Figure 3c). By multiplying the acquired, high- 
resolution data 6 with the reconstruction filter 9, an unfolded image series 11 is 
obtained after Fourier transform along the time axis (Figure 3d). 
[0023] The following describes the data acquisition scheme for a typical 
application of the present invention: 3D single-breath cine SSFP cardiac imaging. In 
addition, the following describes the problem with a straightforward implementation 
and describes the solution to overcome this problem. The solution is part of the 
present invention. 

[0024] Sampling along a sheared grid pattern was implemented into the 
acquisition software of a Philips Gyroscan Intera system (Philips Medical Systems, 
Best, The Netherlands). An elliptical ^-space shutter 12 was applied to reduce the 
number of points to be sampled by approximately 25% by skipping the corners of k- 
space (Figure 4a). The straightforward approach of sampling /c-space is to use a 
linear profile order, which is the default order, but it leads to image artifacts 13 
(Figure 4c). These artifacts originate from the fact that at the end of each kz line 14, 
the kz position is incremented by one, but the ky position now jumps to the starting 
position on the next kz line 15. This will be referred to as a line jump hereafter. 
Usually, when the jump is small, the rewinding gradient is similar in magnitude but 
opposite in polarity to the phase-encoding gradient, so eddy currents generated by 
the former are well compensated by those of the latter. In the case of line jumps, the 
eddy currents are not compensated, resulting in image artifacts. The same also 
applies at the end of each cardiac phase 16. At that point, both the ky and kz 
positions jump back to the starting positions in the next cardiac phase 17. The jump 
is usually small for kz but can be quite large in ky, potentially leading to eddy current 
problems. 

[0025] For a well tuned gradient system, these artifacts should be relatively 
small. Particularly for 2D imaging, these artifacts are almost at a negligible level 19 
(Figure 5a) (see B.A. Jung et. al. Magn. Reson. Med. 2002; 48, p. 921-925). 
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However, these artifacts are drastically amplified in case of accelerated acquisition 
20, particularly at high acceleration factors. This is due to two reasons. Firstly, k- 
space is sampled sparsely in order to accelerate the acquisition. This means that 
larger jumps 21 are involved going from one profile to the next, thereby resulting in 
more residual eddy currents. Secondly, since fewer profiles are acquired at each 
time point or cardiac phase, the line jumps occur more frequently. For example, 
suppose that in real-time 2D imaging, N profiles are acquired. If the usual linear 
profile order is used, a line jump would occur after every N lines. If the acquisition is 
accelerated by, for example, 8-fold, a line jump would occur after every N/8 lines. 
Therefore, the large jumps in /c-space occur more frequently, and the generated 
eddy currents remain significant for the most part of data acquisition, so the acquired 
data are contaminated strongly by image artifacts. 

[0026] To reduce eddy-current related distortions induced by the changing 
phase encoding gradients, the present invention employs a known eddy current 
compensation technique (see B.A. Jung et. al. Magn. Reson. Med. 2002; 48, p. 921- 
925). Previously, this technique was only used for unaccelerated scans, for which 
the eddy current artifacts remain relatively mild for a well tuned hardware system, as 
mentioned. Therefore, it was even mentioned that such a compensation technique 
was unnecessary in the case of 2D imaging (see B.A. Jung et. al. Magn. Reson. 
Med. 2002; 48. p. 921-925). 

[0027] The present invention is further based on the observation that the 
eddy-current artifacts are drastically amplified with accelerated data acquisition. As a 
result, an eddy-current reduction technique, such as the said technique, must be 
used, even in the case of 2D imaging. In the said technique for eddy-current 
reduction, an alternating profile order is used to minimize the large jumps in /c-space. 
For example, in 2D real-time imaging, the profiles are ordered such that /c-space is 
swept along one direction at one time point, and along the opposite direction in the 
next time point (see arrow indicating seep direction in Figures 5c and 5d). This 
profile order significantly reduces the image artifacts both in normal unaccelerated 
acquisition 22 (Figure 5c) and in highly accelerated (e.g. 8-fold) acquisition 23 
(Figure 5d). In cine 3D acquisition, consecutive profiles in the ky-kz plane are 
traversed in a contiguous fashion with the overall direction of sampling being 
reversed for odd cardiac phases (Figure 4b). This alternating profile order scheme 
for subsequent kz profiles and cardiac phases reduces eddy-current related 
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distortions by ensuring maxinrjum connpensation of phase-encode gradient lobes, 
which leads to substantial reduction in image artifacts (Figure 4d) 18. 
[0028] It should be noted that the large jumps in /c-space and the resulting 
eddy-current problems may also affect other imaging applications as well, including 
parallel imaging at high acceleration factors, and reduced field-of-view imaging with 
a highly reduced field of view. Previously, this problem had not been reported in 
those applications because they are generally unable to achieve the high 
acceleration factors used in the present invention. For parallel imaging, this problem 
may begin to arise due to the recent interest in using higher magnetic field strengths. 
At higher field strengths, the achievable acceleration factor is also Increased, so the 
aforementioned problem is expected, and the proposed solution can be applied to 
that application as well. 

Tvpical embodiment: 3D cine cardiac imaging 

[0029] Using 5-fold undersampling (5x k-t BLAST), volumetric data sets 

covering the heart with 20 slices at a spatial resolution of 2x2x5 mm^ were recorded 
with 20 cardiac phases in a single breathhold of 20-22 sec. with the training data 
acquired in a separate 5 sec prescan. When the training data were acquired together 
with the undersampled data in a single scan (referred to as interleaved acquisition 
hereafter), the breathhold duration was prolonged slightly to 25-27 sec. The central, 
densely-sampled region in the ky-kz plane consisted of 49 profiles providing low- 
resolution images without aliasing. A schematic of the sampling pattern applied is 
given in Figure 6. The central portion 15 of the ky-kz space is sampled at the full FOV 
in each cardiac phase whereas the outer portion 16 is under-sampled by a factor of 
five with the pattern being shifted as a function of cardiac phase (Figure 6b). The 
gray dots 24 indicate the sampled points, the dashed dots 25 represent points 
sampled in the previous cardiac phase. Remaining acquisition parameters were as 
follows: FOV: 320x210 mm^ TE/TR: 1.55/3.1 ms, flip angle: 45 deg.. heart phase 
interval: 28-35 ms. Taking the acquisition of low-resolution training data into account, 
the net acceleration factor was 4.3. 

[0030] Navigator-based volume tracking (see e.g. M.V. McConnell et. al. 
Prospective adaptive navigator correction for breath-hold MR coronary angiography. 
Magn. Reson. Med. 1997; 37, p. 148-152) was implemented in order to ensure 
consistency between the low-resolution, fully sampled data and the high-resolution. 
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under-sampled data when the two datasets were acquired in separate acquisitions. 
Using pencil-beam excitation on the diaphragm, the navigator signal corresponding 
to the breath-hold position in the first acquisition was recorded and correlated with a 
navigator acquired prior to the second acquisition. Volume tracking was applied to 
account for possible shifts in breath-hold positions. 

Image reconstruction 

[0031] The data were divided into two sets for reconstruction: a training set 
(Figure 3a) and an under-sampled set (Figure 3c). The training set consisted of the 
densely sampled central region of /c-space only 4; while the under-sampled set 
consisted of all the data sampled along the sheared grid pattern 5. Both data sets 
were inverse Fourier transformed along all axes. The training set yielded a set of 
low-resolution images in x-f space 8. while the under-sampled set yielded a set of 
aliased images in x-f space 6. For each set of aliased voxels in x-f space, the 
aliasing was resolved by the following equation: 



Vector p represents the reconstructed (unaliased) intensities. 6 is a diagonal matrix, 
with the diagonal elements representing the squared expected intensities, as 
determined from the training data. E denotes the encoding matrix, representing the 
process of aliasing. W is the noise variance of the acquisition data. Vector a holds 
the aliased intensity of the voxels. The time-averaged signals of the object (aoc) 
were subtracted prior to reconstruction. After reconstruction, the object 
corresponding to the time-invariant signals (poc) was added back. The square- 
bracketed portion corresponds to the reconstruction filter 9 as shown in Figure 3b. 
The training data were substituted into the reconstructed image for improved data 
consistency, if the training data were acquired in an interleaved fashion. To reduce 
computation, reconstruction was performed separately for each coil element, and the 
reconstructed images were combined afterwards as the root-mean-square of the coil 
images. In general, k-t SENSE, which takes into account the complex sensitivities of 
the receiver coils, may also be used, instead of the coil-by-coil k-t BLAST 
reconstruction described above. 
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In-vivo validation 

[0032] The study population consisted of four healthy volunteers and four 
cardiac patients undergoing cardiac MR imaging for assessment of congenital heart 
disease. In both groups, conventional, multi-slice 2D SSFP acquisitions were 
obtained additionally for comparison, which required 9 sec per slice. Each slice was 
acquired in a separate single breath-hold with imaging parameters comparable to the 
single breath-hold 3D k-t BLAST acquisition. Left ventricular contours were extracted 
using standard post-processing software, and the contour area from the single-slice 
and 3D acquisitions was compared. 

Results 

[0033] 3D data sets were successfully acquired in all subjects. Nearly artifact- 

free 3D volumes were reconstructed, allowing reformatting along arbitrary views. 
Slight image ghosting remained from residual eddy-current effects. In a few subjects 
with large variations of the R-R interval during the breath-hold, early diastolic frames 
exhibited a small degree of blurring. 

[0034] In Figure 7a, a cine 3D single breath-hold data set of central fifteen 
slices at the end systolic phase of a 5x k-t BLAST acquisition obtained in a healthy 
subject is shown. Figures 7b and c depict selected time firames of a mid-ventricular 
short-axis slice (Figure 7b) along with reformatted four-chamber views (Figure 7c). 
The actual temporal resolution was 35 msec. 

[0035] Figure 8 shows example images obtained from one of the patients who 
had repaired Tetralogy of Fallot and pulmonary regurgitation causing right ventricular 
dilation. In all patients, training data were obtained in a separate prescan, lasting 5 
sec, in order to allow shorter breath-hold durations of the high-resolution 3D 
acquisition. The data set was obtained similar to the set of Fig. 7 with 5x k-t BLAST. 
The actual temporal resolution was 28 msec. 

[0036] In in-vivo validation, the mean difference of left ventricular contours 
extracted in three exemplary slices of the 3D data set and the multi-slice 2D 
acquisition was 3.41 ±6.52%. 

[0037] In the present invention it is shown that k-t BLAST can be applied 
successfully to cine 3D acquisition of the heart at 5-fold acceleration. The increased 
scan efficiency allows the acquisition of volumetric data sets of the heart at high 
spatial (2x2x5 mm^) and temporal (28-35 ms) resolutions in a single breath-hold 
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lasting 20-22 sec (for a heart rate of 55-80 beats/min), with the training data being 
obtained separately. The acquisition strategy is well tolerated by healthy subjects 
and patients. 

[0038] The modified profile order according to the invention proved to be 
important even for 2D imaging, particularly when high acceleration factors are used. 
The invention proved effective in reducing eddy-current related image artifacts by 
eliminating the frequent large jumps in ^-space associated with sparse sampling. In 
cine imaging, residual artifacts are considered to be related to the varying number of 
dummy cycles at the end of each cardiac cycle, depending on the R-R interval. 
Further improvements in eddy current reduction can be achieved by using dummy 
cycles after each trigger to provide sufficient time for the eddy currents to decay, or 
to smooth any /c-space jumps. 

[0039] To ensure consistency between breathhold levels when training data 
are acquired in a separate prescan, navigator-based volume tracking is applied. The 
correction factor is set to 0.8, which is considered only a rough estimate of the 
correlation between diaphragm and heart displacement. Inaccuracies, however, do 
not affect reconstruction to a large extend as training data are of very low spatial 
resolution (2x23x1 1 mm^). 

[0040] Fig. 9 shows dtagrammatically a magnetic resonance imaging system 
in which the invention is used. 

[0041] The magnetic resonance imaging system includes a set of main coils 
10 whereby a steady, uniform magnetic field is generated. The main coils are 
constructed, for example in such a manner that they enclose a tunnel-shaped 
examination space. The patient to be examined is slid on a table into this tunnel- 
shaped examination space. The magnetic resonance imaging system also includes a 
number of gradient coils 11, 12 whereby magnetic fields exhibiting spatial variations, 
notably in the form of temporary gradients in individual directions, are generated so 
as to be superposed on the uniform magnetic field. The gradient coils 11, 12 are 
connected to a controllable power supply unit 21. The gradient coils 11, 12 are 
energized by application of an electric current by means of the power supply unit 21. 
The strength, direction and duration of the gradients are controlled by control of the 
power supply unit. The magnetic resonance imaging system also includes 
transmission and receiving coils 13, 15 for generating RF excitation pulses and for 
picking up the magnetic resonance signals, respectively. The transmission coil 13 is 
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preferably constructed as a body coil whereby (a part of) the object to be examined 
can be enclosed. The body coil is usually arranged in the magnetic resonance 
imaging system in such a manner that the patient 30 to be examined, being arranged 
in the magnetic resonance imaging system, is enclosed by the body coil 13. The 
body coil 13 acts as a transmission aerial for the transmission of the RF excitation 
pulses and RF refocusing pulses. Preferably, the body coil 13 involves a spatially 
uniform intensity distribution of the transmitted RF pulses. The receiving coils 15 are 
preferably surface coils 15 which are arranged on or near the body of the patient 30 
to be examined. Such surface coils 15 have a high sensitivity for the reception of 
magnetic resonance signals which is also spatially inhomogeneous. This means that 
individual surface coils 15 are mainly sensitive for magnetic resonance signals 
originating from separate directions, i.e. from separate parts in space of the body of 
the patient to be examined. The coil sensitivity profile represents the spatial 
sensitivity of the set of surface coils. The transmission coils, notably surface coils, 
are connected to a demodulator 24 and the received magnetic resonance signals 
(MS) are demodulated by means of the demodulator 24. The demodulated magnetic 
resonance signals (DMS) are applied to a reconstruction unit. The reconstruction unit 
reconstructs the magnetic resonance image from the demodulated magnetic 
resonance signals (DMS) and optionally on the basis of the coil sensitivity profile of 
the set of surface coils. The coil sensitivity profile has been measured in advance 
and is stored, for example electronically, in a memory unit which is included in the 
reconstruction unit. The reconstruction unit derives one or more image signals from 
the demodulated magnetic resonance signals (DMS), which image signals represent 
one or more, possibly successive magnetic resonance images. This means that the 
signal levels of the image signal of such a magnetic resonance image represent the 
brightness values of the relevant magnetic resonance image. The reconstruction unit 
25 in practice is preferably constructed as a digital image processing unit 25 which is 
programmed so as to reconstruct the magnetic resonance image from the 
demodulated magnetic resonance signals and optionally on the basis of the coil 
sensitivity profile. The digital image processing unit 25 is notably programmed so as 
to execute the reconstruction in conformity with the present invention. The image 
signal from the reconstruction unit is applied to a monitor 26 so that the monitor can 
display the image information of the magnetic resonance image (images). It is also 
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possible to store the image signal in a buffer unit 27 while awaiting further 
processing, for example printing in the form of a hard copy. 

[0042] In order to form a magnetic resonance image or a series of successive 
magnetic resonance images of the patient to be examined, the body of the patient is 
exposed to the magnetic field prevailing in the examination space. The steady, 
uniform magnetic field, i.e. the main field, orients a small excess number of the spins 
in the body of the patient to be examined in the direction of the main field. This 
generates a (small) net macroscopic magnetization in the body. These spins are, for 
example nuclear spins such as of the hydrogen nuclei (protons), but electron spins 
may also be concerned. The magnetization is locally influenced by application of the 
gradient fields. For example, the gradient coils 12 apply a selection gradient in order 
to select a more or less thin slice of the body. Subsequently, the transmission coils 
apply the RF excitation pulse to the examination space in which the part to be 
imaged of the patient to be examined is situated. The RF excitation pulse excites the 
spins in the selected slice, i.e. the net magnetization then performs a precessional 
motion about the direction of the main field. During this operation those spins are 
excited which have a Larmor frequency within the frequency band of the RF 
excitation pulse in the main field. However, it is also very well possible to excite the 
spins in a part of the body which is much larger man such a thin slice; for example, 
the spins can be excited in a three-dimensional part which extends substantially in 
three direcfions in the body. After the RF excitation, the spins slowly return to their 
initial state and the macroscopic magnetization returns to its (thermal) state of 
equilibrium. The relaxing spins then emit magnetic resonance signals. Because of 
the application of a read-out gradient and a phase encoding gradient, the magnetic 
resonance signals have a plurality of frequency components which encode the 
spatial positions in, for example the selected slice. The /c-space is scanned by the 
magnetic resonance signals by application of the read-out gradients and the phase 
encoding gradients. According to the invention, the application of notably the phase 
encoding gradients results in the sub-sampling of the /c-space, relative to a 
predetermined spatial resolution of the magnetic resonance image. For example, a 
number of lines which is too small for the predetermined resolution of the magnetic 
resonance image, for example only half the number of lines, is scanned in the k- 
space. 



